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ABSTRACT 

The light echo systems of historical supernovae in the Milky Way and local group galaxies provide 
an unprecedented opportunity to reveal the effects of asymmetry on observables, particularly optical 
spectra. Scattering dust at different locations on the light echo ellipsoid witnesses the supernova 
from different perspectives and the light consequently scattered towards Earth preserves the shape 
of line profile variations introduced by asymmetries in the supernova photosphere. However, the 
interpretation of supernova light echo spectra to date has not involved a detailed consideration of 
the effects of outburst duration and geometrical scattering modifications due to finite scattering dust 
filament dimension, inclination, and image point-spread function and spectrograph slit width. In 
this paper, we explore the implications of these factors and present a framework for future resolved 
supernova light echo spectra interpretation, and test it against Cas A and SN I987A light echo spectra. 
We conclude that the full modeling of the dimensions and orientation of the scattering dust using the 
observed light echoes at two or more epochs is critical for the correct interpretation of light echo 
spectra. Indeed, without doing so one might falsely conclude that differences exist when none are 
actually present. 

Subject headings: supernovae: individual(Cas A, SN 1987A) — supernovae:general — supernova rem- 
nants — ISM:dust — methods: data analysis — techniques: spectroscopic 



1. INTRODUCTION 
The first light echoes (LEs) were discovered around 



Nova Persei 1901 (|Ritchey|1901a|b[|1902||Kapteynll 902| 
Since then, LEs (whereby we mean a simple scattering 
echo rather than fluorescence or dust re-radiation) have 



been seen in the Galactic Nova Sagittarii 1936 ( |Swope 
1940|) and the er uptive variable V838 Monocerotis (e.g., 



Bond et aL||2003 ) . LEs have also been observed from ex 
tragalactic supernovae (SNe), with SN 1987A (?) being 
the most famous case (Crotts||1988| |Suntzeff et al.||1988| 
Bond et al |1990l|Xu et al|l()95NSugerman et al j2M5a|b| 
Newman fc Rest||2006D, but alsoTnclud ing SNe 1991T 
([Schm idt et al. 1994; Spark s et aL"1|1999[) , 1 993J ([Suge~ 
man & CrotTsporjgfjLiu et al.|20U3t"T9" g5E ( iQuinneFar 
2006|), 1998 bu (Cappell aro et al.||2001|), 2002hh flMe jHe 
eTar||2006[ |Wel ch||2007|), 20 U3gd jSugerman|l 2005; Van 
Dyk et ai.||2' 



6), and 2U06X flWang ; et al.||2008| ). 

The suggestion that historical SJNe mi ght be studied 
by their scattered LEs was first made by Zwicky (19401 



and attempted by van den Bergh 



LEs of centuries-old SINe were d iscovere d in t he Large 



(1965 1966). TheTirst 



Magellanic Cloud (LMC) ( |Rest et al.| |2005b[ ). They 
found three LE complexes associated with three small 
and therefore relatively young SN remnants (SNRs) in 
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the LMC which were subsequently dated as being be- 
tween 400-1000 years old using the LMC di stance and 
LE fea ture positions and apparent motions (Rest et al. 



2005b). These findings provided the extraordinary op 



portunity to study the spectrum of the SN light that 
reached Earth hundreds of years ago and had never been 
recorded visually or by modern scientific instrumenta- 
tion. A spectrum of one of the LEs associated with 
SNR 0509-675 revealed that t he reflected LE light came 
from a high-l uminosity SN la ([Rest et al. | 2008a[), similar 
to SN 1991T ( |Filij22enko_et al.l 1 992 [ | P hil'ips et al.|l992 | 
and SN 1999aa ( |Uaravini et al.||2U04[ ) — the first time 
that an ancient SJNe was classified using its LE spec- 
trum. Analysis o f ASC A X-ray data of SNR 0509-675 



led Hughes et al. ( 1995 ) to suggest that this SNR likely 
came from a S IM la. Recent analysi s of its Chandra X- 
ray spectra by Badenes et al. ( 2008 ) also supported the 
classification of this object as a high- luminosity SN la. 

The inte rpretation of LEs, h owever, is fraught with 
subtleties. |Krause et al. (2005), for example, identified 
moving Cas A features (called infrared echoes" or "IR 
echoes") using mid-IR imaging from the Spitzer Space 
Telescope. IR echoes are the result of dust absorbing the 
outburst light, warming and re-radiating a t longer wave- 
length s. The main scientific conclusion of |Krause et al.| 
(2005) was that most of the IR echoes were caused by 
a series of recent X-ray outbursts from the compact ob- 
ject in the Cas A SNR based on their apparent motions. 
However, the analysis was flawed because it did not ac- 
count for the fact that the apparent motion strongly de- 
pends on the inclinatio n of the scattering dust fila ment 
( |Dwek fc Arendt||2008| ?). |Dwek fc Arendt| ( |2008[ ) also 
showed that X-ray flares are not energetic enough to be 
the source of these IR echoes, but that instead the LEs 
must have been generated by an intense and short burst 
of EUV-UV radiation associated with the break out of 
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Rest et al. 



the shock through the surface of the Cas A progenitor 
star. Therefore, the most likely source of all LEs asso- 
ciated with Cas A, both infrared and scattered, is the 
Cas A SN explosion itself. The need to properly analyze 
the growing collection of LE spectra is what motivated 
the present work. 
The hrst scattered LEs of Galactic SNe associated with 



Tycho ' s SN and the Cas A SN were disco vered by|Rest 
et al. (2007, 2008b). Contemporaneously, Krause et al. 
Q20U8p obtained a spectrum of a scattered optical LE 
spatially coincident with one of the Cas A IR echoes, 
and identified the Cas A SN to be of Type lib. Because 



Rest et al. (2008b| discovered several LEs with multiple 
position angles relative to the SNRs, there is a new op- 
portunity for observational SN research — the ability to 
measure the spectrum of the same SN from several differ- 
ent directions. At any given instant after the SN light has 
reached an observer by the direct path, an ellipsoidal sur- 
face exists where scattered light from the SN may reach 
the observer, whose arrival to the observer is delayed by 
the additional path length. Interstellar dust concentra- 
tions must lie at one or more locations on the ellipsoid 
for the scattered LE to be detectable. Since each LE is at 
a unique position on the ellipsoid, each LE has a unique 
line of sight to the SN. Therefore, the scattered spectral 
light can provide spectral information on the asymmetry 
of the SN photosphere. In the limiting case of dust on 
the ellipsoid opposite the observer, the spectrum of the 
scattered LE carries the signature of conditions on the 
SN hemisphere usually hidden from the observer. As the 
time since the explosion increases, the scattering ellip- 
soid surface expands and may intersect additional dust 
concentrations. We present a framework on how to in- 
terpret observed LE spectra depending on the scattering 
dust properties, seeing, and spectrograph slit width and 
position. 

In this paper, we examine the pitfalls of previous analy- 
ses, suggest improvements for future studies, and present 
the application of our methods to data. In Section [2] we 
describe how the properties of the scattering dust can 
cause significant differences to observed LE spectra, and 
show a specific example of how this has been overlooked 
in the past. In Section [3] we quantify how observed LE 
spectra depend on potential differences in dust character- 
istics (thickness and inclination), as well as observational 
characteristics (seeing and spectroscopic slit size and ori- 
entation) . We then apply this framework to observations 
of Cas A (Section [§ and SN 1987A (Section ||. In Sec- 
tion [6j we discuss now differences in the dust and ob- 
serving conditions can effect an analysis of LE spectra 
and suggest future applications of LE observations. We 
conclude in Section 

2. A CASE STUDY: THE KRAUSE ET AL. (2008) CAS A 
LIGHT ECHO SPECTRUM 

A common assumption is that an observed LE spec- 
trum is equivalent to the light-curve weight ed integra- 
tion of the spectra at indi vidual epochs (e.g. jRest et al. 
2008a Krause et al. 12008 ). However, Figure [l] illustrates 
how slit width, dust filament width and inclination signif- 
icantly influence the observed LE spectrum. We consider 
a LE from a 300-year-old event 10,000 light years away 
that lasts 150 days (yellow shaded area), with a rising 
and declining arm of 30 and 120 days, respectively. The 



peak of the event is indicated with the red line. We place 
the origin of our coordinate system p, z at the source 
event. The z-axis is the distance of the dust to the source 
event along the line of sight. The p-axis is in the plane 
of the sky from the source event to the LE position. The 
quadratic relation between the parameters z, p, and the 
time since explosion t is described by the well-known LE 
equation ( Couderc||1939 ) (see Section [A] for details). 

The scattering dust filament (brown shaded area) is lo- 
cated at (p, z) = (300, 0) ly with an inclination of a = 0° 
where a is the a ngle with respect to the plane of the 
sky (see Equation [A3| . The dust width is 0.02 ly, 0.2 ly, 
and 0.6 ly from left to right, respectively. The right- 
most panel of Figure [l] shows the flux of the LE that we 
would observe with these dust filaments versus p. We de- 
note these as the LE profiles, which essentially are cuts 
through the LE along the axis pointing toward the source 
event. In order to be able to directly compare the LE pro- 
files to the dust and LE pulse geometry, we show the flux 
and p as the x- and y-axis, respectively. Note that for 
very thin dust filaments like the one with ad = 0.02 ly, 
the LE profile is the the projected light curve. 

A spectroscopic slit with a width of 1" is indicated with 
the gray shaded area. In the leftmost panel, the dust 
filament has a dust width of 0.02 ly (0.006 pc). Only 
light that is in the intersection of the event pulse, dust, 
and slit (indicated with the thick black rectangle) con- 
tributes to the observed spectra. In this case, only 20% 
of the event (30 out of 150 days, indicated by the dashed 
blue lines) contribute to the spectrum. The second panel 
from the left shows the same example, but this time the 
dust filament has a width of 0.2 ly (0.06 pc), and ~51% 
of the event pulse (77 out of 150 days) contributes to 
the spectrum. Only if the dust width is at least 0.6 ly 
(0.18 pc) thick (third panel from the left), the full event 
contributes to the spectrum. 

This clearly demonstrates that the observed LE spec- 
trum is not necessarily the full light-curve weighted in- 
tegrated spectrum when the scattering dust filament is 
sufficiently thin. The effect of a thin filament of dust on 
the interpretation of a LE spectrum is best appreciated 
by considering the difference in persistence of spectral 
features at early and late time. Naively, one might ex- 
pect that the faintness of a SN at late times would not 
result in spectral features from those phases contributing 
to the integrated spectrum in a significant way. However, 
in the late nebular phases, most of the flux is typically 
concentrated in very few lines. The top panel of Fig- 
ure [2] portrays the spectra of SN 1993J, the very well- 
observ ed and prototypical exa mple of the SN lib class 



(e.g., ^Richmond et al. 1994 



__?;?), at 6, 39, 88, and 
205 days after peak brightness from the bottom to the 
top, respectively (???). During the late phases, a sig- 
nificant fraction of flux is radiated in the [O I] AA6300, 
6364 doublet, which completely dominates the spectra 
for times later than 200 days. The bottom panel shows 
the light-curve weighted integrated spectra of SN 1993J 
with different integration limits. The integration limits 
of these spectra are from 20 days before peak brightness 
to -Pmax, where P max ranges from 5 to 300 days after 
peak brightness from the bottom to the top. The most 
significant difference between the integrated spectra with 
early and late integration limits is the emergence of the 
[O I] doublet. We use the SN 1993 J lightcurve from ? as 
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Fig. 1. — Illustration of how the dust filament width, a^, influences the observed LE spectrum. The yellow shaded area indicates a 
300-year-old SN event that lasts 150 days, with a rising and declining arm of 30 and 120 days, respectively. The peak is indicated with a 
red line. The scattering dust filament (light grey shaded region) is located in the plane of the sky crossing the event position (z = 0.0 ly), 
with an inclination of a = 0°. The dust width is 0.02 ly, 0.2 ly, and 0.6 ly from left to right, respectively, and the left panel shows the 
respective LE profiles. A spectroscopic slit with a width of 1" is indicated with the dark gray shaded area. Only light that is within the 
intersection of the event pulse, dust, and slit (indicated with the thick black rectangle) contributes to the observed spectrum. The blue 
dashed lines indicate which part of the light pulse is probed by the spectrum. 



weights for the integration. 

SNe lib are not unique in this regard. The nebular 
spectra of all SNe are dominated by strong emission lines. 
Specifically, SNe la have strong [Fe II] and [Fe III] lines, 
SNe Ib/c have strong [O I], [Mg I], Ca II, and [Ca II] 
lines, and SNe II h ave strong hydrogen Balmer lines. See 
Filippenko] ( 1997 1 for a review of SN spectroscopy. 

T he perils ot not properly accounting for differences in 
dust filament thickness are clear — the strength of fea- 
tures in observed LE spectra can be significantly altered 
by the width of the scattering dust on the sky. 

In Figure [2[ we di splay the LE spectrum of Cas A by 
Krause et al. (20081 in red. This spectrum has a weak 
[U 1J doublet, consistent with the integrated spectra of 
SN 1993J, but only if integrated until ^80 days after peak 
brightness. The full integrated spectra of SN 1993J have 
a much stronger [O I] feature. If the LE spectrum is the 
full light-curve weighted integrated spectrum of Cas A, 
then Cas A and SN 1993J have a significant difference 
at wavelengths near the [O I] feature. The alternative 
is that the LE spectrum does not contain light from the 
late-time portion of the SN light curve, i.e., the weighting 

function is truncated at late times. 

In their co mparison to the observed Cas A LE, |Krause| 
et al. ( |2008[ ) also used an integrated SN 1993J spectrum 
- but only integrated to 83 days after peak brightness. 
This seemingly prudent choice was not explained in the 
text, but presumably they found this integration time 
to yield the best comparison. However, this analysis 
leaves an open question: Was the Cas A SN different 
from SN 1993J in this regard or are the differences due 
to the properties of the scattering dust? In the follow- 
ing sections, we attempt to answer this question and its 
more general counterpart for all LEs. 



3. LIGHT ECHO PROFILE MODELING 



In Appendix [XI we derive that the observed LE profile 
shape, P, strongly depends on the dust thickness, dust 
inclination, and seeing. We also show that this may have 
a significant impact on the observed LE spectra (see Ap- 
pendix [B^ and O). Here we investigate the magnitude of 
each of these effects on observed LE spectra. 

For this analysis, we examine a toy model of a LE that, 
despite its simplicity in some respects, appears to be rep- 
resentative of actual LE situations. The toy model is a 
LE from a 300-year-old event 10,000 ly from the observer. 
For the event light curve, we choose a "sail" shape, with 
the peak at phase 0.0 years and then decreasing to zero 
flux at a phase of 0.3 years and —0.1 years (see the dot- 
ted curves in the right-hand panel of Figures [3] through 
[7 for a representation of this light curve), sharing the 
characteristic short rise and long decline of a SN light 
curve. The scattering dust filament intercepts the LE 
paraboloid at (po,z ) = (300,0) ly. For the dust model, 
we use a dust filament that is locally planar and perpen- 
dicular to the (p, z) plane, with a boxcar width profile. In 
the following subsections, we vary the properties of a sin- 
gle parameter while holding the others constant. We are 
therefore able to determine the effects of each parameter 
on the observed LE and its spectrum. 

3.1. Dust Thickness 

As described above, the dust thickness can significantly 
alter the LE spectrum. Using our toy model, we can 
quantify this effect. 

Figure [3] presents LE profiles, window functions, and 
effective light curves for the toy model with different dust 
widths. The left column shows the LE profiles for each 
dust thickness. The ir-axis is Ap = p — p , where (po, ^o) 
is the position where the dust intersects with the LE 
paraboloid at phase 0.0 years of the event. We choose a 
dust filament inclination of a = 45° , since it is perpendic- 
ular to the LE paraboloid and thus most representative 
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Fig. 2.— Top panel: Spectra of SN 1993J at +6, +39, +88, 
and +205 days relative to peak brightness from the bottom to the 
top, respectively. During the late phases (> 40 days after peak), 
a significant fraction of flux is emitted in the [O I] AA6300, 6364 
doublet, which dominates the spectra after ~200 days. The [O I] 
doublet is indicated by the yellow-shaded region. Middle panel: 
light-curve weighted integrated spectra of SN 1993J. The integra- 
tion limits of these spectra are from —20 days relative to peak 
brightness to P max , where P max ranges from +5 to +300 days 
from the bottom to the top. The most significant difference be- 
tween the integrated spectra with early and late integration limits 
is the emergence of the [O I] doublet. The Cas A LE spectrum of 
|Krause et al.|p008| l is shown in red, which does not show a promi- 
nent [U 1J doublet feature. Bottom panel: light-curve weighted 
integrated spectra of SN 1993J in the wavelength region near the 
[O I] doublet and Ho features. The integrated spectra have in- 
tegration limits of +80 days and +300 day s are compared to the 
Cas A LE spectrum of |Krause et al.] l |2008| red curve). The inte- 
grated SN 1993J spectrum witn tne +SU day integration limit is a 
better fit to the [O I] doublet in the observed Cas A spectrum. 



of typical inclinations. Centering a 1.0" spectroscopic 
slit on the peak of the LE profile, we calculate the as- 
sociated window function (middle column) and effective 
light curves (right column), which are simply the light 
curve convolved with the window function. The different 
rows show different dust widths. 

There are two distinct regimes where a modest change 
in the dust width has little effect on the effective light 
curve or integrated spectrum. The first one is shown 
in the top panels, where the dust width is so thin that 
the observed LE profile is effectively the projected initial 
light curve (corresponding to the left panel of Figure [l]) . 
The other limit is a very thick dust sheet (corresponding 
to the right panel of Figure [T]). 

In the former case (top row of Figure , changing the 



1.0 
0.8 
<0.6 

3 

; 0.4 

0.2 

0.0 

1.0 

0.8 

<0.6 
i 

; 0.4 
0.2 
0.0 
1.0 
0.8 

<0.6 

3 

: 0.4 
0.2 
0.0 



light ech 
profile 






effective — =0.001 
light \ — ^=0.01 
light 
curve 





-6 -4-2 2 
Ap in arcsec 



-100 100 
phase in days 



-50 50 100 
phase in days 



Fig. 3. — Modeled LE profiles (left column), window functions 
(middle column), and effective light curves (right column) for our 
toy model as a function of dust width. All models in all rows have 
the same parameters except for the dust width; the separation into 
three rows is only for clarity. Ap is with respect to the position 
where the dust intersects the light pulse. The toy model is a 300- 
year-old event with a "sail" shape as described in the text. The 
scattering dust filament intersects the event pulse paraboloid at 
z = 0.0, and has an inclination of a = 45°. 

dust width from a c [ = 0.001 ly to a 4 — 0.01 ly only 
marginally modifies the effective light curve. Centering 
a spectroscopic slit on the LE profile results in a box-car 
window function, and the effective light curve is trun- 
cated to the region near maximum brightness. In the 
latter case (bottom row of Figure 3]), the dust width is 
large enough that the LE profile is dominated by the dust 
width. The LE profile has a flat top, and only the edges 
are impacted by the original light curve profile. Center- 
ing a spectroscopic slit on the LE profile results in very 
wide window functions, and the effective light curve is 
essentially the same as original light curve (see the right 
panel of Figure [I] for an illustration). 

As we show in Sections H] and [5j our observations of 
LEs have been a mixture of both regimes. The middle 
row of Figure [3] shows the light curves, window functions, 
and effective light curves for dust widths of ad — 0.1 ly 
(0.03 pc) and a d = 0.2 ly (0.06 pc). The effective light 
curve is still heavily weighted to the time near maxi- 
mum^] but with increasing dust width the effective light 
curve slowly transitions to the full light curve. 

3.2. Dust Inclination 

The dust inclination has a large effect on the LE pro- 
file. This is simply a projection effect; as the dust sheet 
rotates, the projection of the light curve on the sky can 
change significantly. Turning to our toy model, we quan- 
tify this effect. 

We fix the dust width to o& = 0.1 ly. Figure |4l shows 
modeled LE profiles, window functions, and effective 

8 In this example, with increasing dust width, the peak of the LE 
profile moves away from Ap towards the negative direction. This 
is caused by using a box-car profile for the dust width profile. This 
effect is minimized if the dust width profile is peaked (e.g., if it is 
a Gaussian). 
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Fig. 4. — Same as Figure[3] but for varying dust filament inclina- 
tions, ct, and constant dust width. The width of the dust filament 
is (Tjj = 0.1 ly. Note that the LE profile gets "stretched" when the 
inclination of the dust filament approaches the tangential angle of 
the LE ellipsoid, which is a = —45° for z = ly. 

light curves for several examples that only differ by the 
inclination of the scattering dust filament. As expected, 
we find that the filament inclination has a profound im- 
pact on the LE profile width. For inclinations close to 90° 
(see a = 70° in the top row of FigureEl), the dust filament 
is aligned along the line of sight, while the light curve 
projected into p space is compressed. Consequently, the 
window function has large wings, and the effective light 
curve is very similar to the full light curve. However, for 
inclination angles corresponding to dust sheets close to 
the tangent plane of the LE paraboloid (for z Q = 0.0 ly is 
—45°), the LE profile subtends a larger angle on the sky 
(see a = —30° in bottom row of Figure El). As a result, 
a spectroscopic slit will cover a very smalT portion of the 
light curve, and the window function is approximately a 
box car with a small width. 

It is clear that the determining the dust inclination is 
very important, but fortunately it is straightforward to 
do so from the observations in cases where the location of 
the SN and date of outburst are known. In Appendix [E] 
and in more detail in Rest et al. (in prep.), we derive the 
dust inclination from the apparent motion of the LE in 
difference images, and discuss its sources of uncertainty. 
We discuss the impact the dust inclination has on the 
possibility of sampling the spectra within different age 
intervals for the same SN, as well as the possibility of con- 



straining the SN light curve shape in Section [672] and [673 
respectively. 

3.3. Point-Spread Function Size 

Observations of resolved sources must always account 
for the point-spread function (PSF). LEs are no excep- 
tion since degradation of an image by a large PSF will 
not only smear the LE across the sky, but will also smear 
the LE in time. Again, we turn to our toy model, and 
vary the PSF to determine its effect on the interpretation 
of LEs. 

Figure [5] shows our toy model with dust inclination 
of a = 45°. We examine a Gaussian PSF with four 
different full- widths at half-maximum (FWHMs): 0.0" 
(corresponding to a 5 function PSF), 0.4", 1", and 2". 
Each row of Figure[5]represents a different dust thickness, 
with the top, middle, and bottom rows corresponding 
to <7<2 = 0.001, 0.1, and 0.3 ly, respectively. For small 
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Fig. 5. — Same model as in Figure [3] but with varying the dust 
filament width and seeing. The dust width presented are <r d = 
0.001 ly (top row), cr^ = 0.1 ly (middle row), and = 0.3 ly 
(bottom row). Seeing of 0.0", 0.4", 1.0", and 2.0" is indicated 
with red, blue, yellow, and green lines, respectively. 

dust widths (< 0.1 ly), and thus thin LE profiles, the 
effective light curve depends significantly on the seeing. 
However, for larger dust widths, the seeing has very little 
effect on the effective light curve. This insight provides 
guidance for prioritizing targets for spectroscopic follow- 
up observations based on the conditions at the telescope. 

3.4. Slit Offset 

Thus far, we have assumed for our toy model that the 
spectroscopic slit was centered on the peak of the LE 
profile. However, LEs are usually both very faint and ex- 
tended. Additionally, their apparent motion on the sky 
makes perfect alignment of a spectroscopic slit unlikely - 
especially when there is a significant delay between the 
time of mask preparation and observations. Figure [6] 
shows the effect of having the slit offset from the LE 
profile peak. As before, we choose a representative dust 
inclination of a = 45°, and the three rows, from top-to- 
bottom, correspond to dust widths of ad = 0.001, 0.1, 
and 0.3 ly, respectively. The red, blue, and green shaded 
areas indicate 1" slits with 0.0", —1.1", and —2.2" offsets 
from the peak, respectively. Clearly a 1" offset has a sig- 
nificant impact on the effective light curve for all but the 
thickest dust widths (but even in that case, it could still 
result in a significant difference if the spectroscopic fea- 
tures in the outburst spectra are changing rapidly with 
time) . Figure \f\ shows the same scenario as Figure km but 
with the addition of a 1" PSF to the model. 

Figure [8] illustrates the application of this model to a 
"real world" set of conditions. We define a coordinate 
system (Ap,u), where as before Ap = p — p$, and u is 
orthogonal to the p-axis through Ap — 0.0 in the plane 
of the sky. For a dust filament that is locally planar and 
perpendicular to the (p, z) plane, i.e. parallel to u, the 
LE is aligned with the u axis. However, if it is tilted 
with respect to u by an angle of 5\, then the LE will 
also be tilted by the same angle (see gray shaded area 
in illustration). We consider a slit (red box in Figure [8| 
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Fig. 6. — Same as FigurcR?] but with slits that are offset by 0.0", 
-1.1", and -2.2" from the peak of the LE profile. The 1" slits are 
indicated with the shaded areas in the left column. 




-2 2 
Ap in arcsec 



-100 100 200 
phase in days 



50 50 100 
phase in days 



Fig. 7. — Same as Figure H but in 1 



seeing. 



that is misaligned with the LE by an angle &2, and also 
offset along the p axis by As. 

In this situation, we can recover the proper LE window 
functions. To do so, we can split the slit into multiple 
subslits, (in this example, we show this with the horizon- 
tal dotted lines). For each of these subslits, we determine 
the LE profile, calculate the offset of the subslit to the 
LE using 82 and As as described in Appendix [Cj and de- 
termine the window function for each subslit. The total 
window function for the full slit is then the flux-weighted 
average of the subslit window functions. 

4. CAS A LIGHT ECHOES 

4.1. Observations 

As a case study, we investigate LEs discovered in a 
multi-season campaign beginning in October 2006 on the 
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Fig. 8. — Toy model of a LE and spectroscopic slit configuration 
where the slit is misaligned both in angle and offset relative to 
the LE. The LE is indicated with the shaded area. We define a 
coordinate system (p, u) in the plane of the sky centered at the 
LE, where the p axis is in the direction to the source event, and 
u is perpendicular to p, with the positive direction along positive 
position angle of LE-SNR. The LE is along the u axis if the dust 
filament is perpendicular to the (p, 2) plane. In this example he 
angle 81 between the LE and the n-axis is 18°, and the angle 62 
between slit and LE is 22°. Six subslits of width 0.8" are indicated 
with the blue dotted lines. The offset As is 1.0". 



Mayall 4 m telesco pe at Kitt Peak National Observatory 
Rest et al.||2008b[ ). The MOSAIC imager on the Mayall 
m telescope, which operates at the f/3.1 prime focus at 
an effective focal ratio of f/2.9, was used with the Bern- 
stein "VR" broad filter (A c = 595 nm, 5X = 212 nm; 
NOAO Code kl040). The imaging data was kernel- and 
flux-matched, al igned, subtracted , and masked using the 
SMSN pipeline ([Rest et al.|2005a|[Garg et al.|2007lMrF] 
naitis et al. 2007 T- We investigate the three C as A SIN 
light echoes LE21T6, LE2523, and LE3 923, for which we 
presented spectra in Rest et al. (2010). Figure [9] shows 
the difference images and the slit position for these three 
LEs, and Table[T]lists the positions, lengths, and position 
angles for each slit . LE2116 was previously reported by 
Rest et al.| (|2008b), whereas LE2521 and LE3923 were 
discovered on 2009 September 14 and 16, respectively. 
For the analysis of LE2116, we use images from 2009 
September 14, while we use the discovery images for 
LE2521 and LE3923. 



4.2. Cas A Light Echo Template Spectra 

There are several LEs associated with the Cas A SN 
Here we examine three LEs for whic h we have spectra 
Their profiles are shown in Figure 10 To fit the LE pro- 
files, we numerically solve the model defined by Equa- 
tion |A14| for the observed data. To do this, we a ssume 



that the exp losion occurred in the year 1681 ± 19 (Fesen 



et al. 2006), determine the scatt ering dust coordinates 
Po and zq by using Equation ~Kl with the measured an- 
gular separation between the LE and the SNR, and we 
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RA 


Dec 


PA SNR-LE 


Seeing"" 1 


Width a 


Length 3, 


p A a 


LE 


Telescope 


(J2000) 


(J2000) 


(°) 


(") 


(") 


(") 


(°) 


LE2116 


Keck 


23:02:27.10 


+56:54:23.4 


237.83 


0.81 


1.5 


2.69 


0.0 


LE2521 


Keck 


23:12:03.86 


+59:34:59.3 


299.11 


0.79 


1.5 


4.32 


41.0 


LE3923 


Keck 


00:21:04.97 


+61:15:37.3 


65.08 


0.89 


1.5 


7.72 


0.0 



Parameters for spectroscopic slit 



|LE2116| i 

15" 









Fig. 9. — Difference image cutouts for LE2116, LE2521, and 
LE3923 from left to right, respectively. For all image stamps north 
is up and east is left. Excess flux in the first and second epoch is 
white and black, respectively. The position and size of the spec- 
troscopic slit is indicated by the white rectangle. 
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Fig. 10.— LE profiles for LE2116, LE2521, and LE3923 from top 
to bottom, respectively. The model fit is indicated with the black 
lines. 

determine the dust inclination from the apparent mo- 
tion of the LEs (this process is defined in Appendix [e] 
and described in more detail by Rest et al., in prep.jT^ 
The remaining free model parameters are the LE profile 
peak height and position and the dust filament width, 
which has physical implications. For the three LEs pre- 
sented above, we display their best fit model profiles in 
Figure [10} 

To numerically fit the model parameters, we determine 
the best fitting profile by calculating the reduced % 2 for a 
3-dimcnsional grid of the three fit parameters. For each 
parameter, we marginalize over the other parameters to 
obtain the best-fitting value. The initial step sizes of the 
parameter grid are often too large for a proper marginal- 
ization. To ensure a proper determination of the param- 
eters, we iterate this process three times adjusting the 
grid spacing, guided by the fitted parameters and their 
uncertainties from the previous iteration. In all but the 
lowest signal-to-noise ratio cases, the probability dens ity 
functions have a well-behaved Gaussian shape. Figure 
shows as examples the marginal probability density func- 
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Fig. 11. — Probability density functions (PDFs) for the peak 
position (left) and dust width (right) of LE2521 ata = 0.5". The 
most likely value is indicated with a dashed line, and the yellow 
shaded region indicates the 1-<t interval. 



tions (PDFs) for the peak position and dust width of the 
LE2521 LE profile at u = 0.5". 

Table 2] shows the input and fitted parameters of these 
three LE profiles. These profiles are for subslits with 
u = 0.5", i.e., close to the center of the spectroscopic slit. 
The definition of the {p, u) coordinate system and how it 
relates to the LE and the subslits is explained in detail 
in Section [374], Appendix [C] and Figure [§] LE2116 and 
LE2521 have profiles with a similar width. The profile of 
LE3923, however, has a significantly larger width. Note 
that the model fits the data very well. 

Figure [12| shows the fitted LE peak position p , dust 
width ad, and peak surface brightness SB for all subslits 
of each of the three LEs. Panels A-C show po for the 
three LEs, with the spectroscopic slit overplotted (black 
solid lines) . The dashed line indicates a straight line fit to 
the subslits that are within the spectroscopic slit. Even 
though there is some structure in the po versus u rela- 
tion, locally the LE position can be approximated very 
well with a straight line. Panel D shows the dust width 
o-rf. The dust width of LE3923 is nearly constant over 
20" and is an order of magnitude larger than the dust 
width of the other two LEs. The dust width for LE2116 
decreases and then increases over only ~2". The width 
decreases to ad = 0.03 ly for u > 0.5, making this portion 
of the LE2116 filament the thinnest yet observed. We dis- 
cuss potential scientific opp ortun ities for a LE with such 
a thin dust width in Section ROl Panel E shows the sur- 
face brightness. LE2521 has a surface brightness brighter 
than 21 mag/arcsec -2 and is the brightest Cas A LE yet 
reported. However, it is spatially small. On the other 
hand, LE3923 is relatively faint, but its large size makes 
spectroscopy productive. 
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TABLE 2 
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Seeing 
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5i 


s 2 
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Pa 
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SB 


Cs,eff 


LE 


UT Date 


(") 


n 


(°) 


(°) 


(") 


(iy) 


(iy) 


(iy) 


(mag) 


(") 


LE2116 


20090914 


0.97 


9(5) 


7.80 


-39.97 


-0.42 


627.9561(18) 


437.11 


0.031(10) 


22.82(02) 


1.77 


LE2521 


20090914 


1.13 


54(5) 


-10.98 


-0.91 


-0.01 


317.9613(23) 


-9.89 


0.093(04) 


21.83(02) 


1.53 


LE3923 


20090916 


1.65 


7(5) 


10.44 


-35.36 


-0.35 


1203.8904(76) 


2045.38 


0.614(53) 


23.35(04) 


1.65 



Note. 



All the values arc for the subslits with u — 0.5. 
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Fig. 12. — Panel A, B, and C show the positions of the LE profile 
peaks for the different subslits for LE2116, LE2521, and LE3923, 
respectively. Apo(u) is defined as po(u) — po(0), and u is orthogonal 
to the p-axis through po(0). Red, blue, and cyan symbols indicate 
LE2116, LE2521, and LE3923, respectively. Panel D shows the 
fitted dust widths, cr^, and Panel E shows the surface brightnesses. 



For the spectra presented by Rest et al. (2010), the 
profile fits had not been completed when slit masks were 
designed. As a result, the alignment between the LEs 
and the slit were not optimal. In particular, there is a 
large misalignment in angle for LE3923; however, that 
LE has a larg e profile (FWHM » 10"; see bottom panel 
of Figure 10 1 and therefore an offset at the 1-2" scale 
does not decrease the total amount of flux by a signifi- 
cant amount. However, it does affect the observed spec- 
tra since the window function and effective light curves 
strongly depend on these offsets. 

For each LE, we calculate the window function by fit- 
ting the observed LE profile with the numerical model as 
described above. The spectroscopic slit is applied to the 
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Fig. 13. — Top: window function for the different LEs. Middle: 
effective light curves of SN 1993J for the LEs. The unmodified 
light curve of SN 1993J is shown in black. Bottom: integrated 
SN 1993J spectra, where the integration is weighted by the respec- 
tive effective light curve. Red, cyan, and blue curves represent 
LE2116, LE2521, and LE3923, respectively. The black spectrum 
is weighted by the original, unmodified SN 1993J light curve. The 
differences in the spectra are mainly due to the strong [O I] AA6300, 
6363, [Ca II] AA7291, 7324, and Ca II NIR triplet features in the 
late-phase spectra. 



fitted LE profile using Equations |Blf|B4| as described in 
Appendix [Bl It is important to note that now the see- 
ing from the spectroscopic observation is used instead of 
the seeing from the image used to determine the LE pro- 
file. For each slit, the weig hted aver age of all subslits is 
calculated using Equatio ns |C1| - C4 in Appendix [Cj 



The top panel of Figure [l3| shows the window functions 
for the three LEs for whichwe have spectra. The middle 
panel shows the SN 1993J light curve and the effective 
light curves for each LE (using the LE profiles and as- 
suming that their light curves are the same as SN 1993 J). 

Even though LE2521 and LE2116 have comparable 
dust widths, the dust inclination for LE2521 is signifi- 
cantly larger, resulting in a wider window function than 
LE2116. This also causes the projected SN light-curve 
shape to be "squashed." LE3923 has a similar inclina- 
tion to that of LE2116, but its scattering dust is thicker, 
also leading to a wider window function. Note that 
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the three window functions, and therefore also the ef- 



fective light curves (see middle panel in Figure 13), have 
significantly dropped off after 100 days past maximum. 
However, the decline time of SN is significantly longer, 
and consequently the contribution of spectra with phases 
> 100 days is cut off - a fact illustrated in the bottom 
panel of Figure[l3j which shows the integrated spectra for 
SN 1993 J. The late type spectra of SN lib like SN 1993 J 
are dominated by [O I] AA6300, 6363, [Ca II] AA7291, 
7324, and the Ca II NIR triplet. Note that the black 
spectrum, which is weighted by the original, un-modified 
light curve, is much stronger at the position of these lines 
compared to the other spectra, which are integrated us- 
ing the effective light curves. We conclude that a per- 
fectly symmetrical SN observed at different position an- 
gles might be erroneously interpreted as asymmetrical if 
scattering dust properties, seeing, and slit width are not 
taken into account. 

4.3. Light Echo Spectra: Systematics 

We have shown in Scction|3]with simulated light curves 
and dust filaments how the LE spectra are influenced by 
the dust properties and observing conditions, and in Sec- 
tion |4.2| we have applied this to the real- world example 
of the Cas A LEs. In this section, we investigate how 
changing these properties impacts the LE spectra. We 
start with the fit to the Cas A LE profiles of LE2521 
described in Section |4.2| and create a window function, 
effective light curve, and LE spectrum using using the 
spectroscopic and photometric library of SN 1993J as a 
template. We then vary dust width, dust inclination, 
and seei ng w hile keeping all other parameters fixed. 

Figure 14 shows the effective light curves (left panels) 
and integrated SN 1993 J spectra (right panels). In the 
top panel, we vary the dust width, We find that, 
consistent with the simulations in Section [3] the nebular 
[O I] and [Ca II] lines in the integrated spectra have a 
strong dependence on the dust width. For ad % 0.1 ly, 
the [O I] and [Ca II] emission lines are weak, whereas for 
larger values the integrated spectrum quickly approaches 
the limit of the integrated spectrum without any window 
function applied with strong [O I] and [Ca II] emission 
lines. We find similar results for the dust inclination: the 
difference between the integrated spectra is small and the 
[O I] and [Ca II] lines are weak for a < 60°, but then 
rapidly approaches the integrated spectrum without any 
window function applied for a close to 90° . We find that 
the seeing (bottom panel) only has a very strong impact 
on the observed LE spectru m if the dus t width is very 
small as discussed in Section [6~2l and [Ol 

We also investigate how uncertainties in the dust in- 
clination impact the derived integrated spectra template. 
To do this, we fixed the dust inclination from 0° to 80°for 
LE2521, and fit the LE profile. Figure [15] shows the fitted 
dust width, a,i, with respect to the input dust inclination, 
a. Three examples of t he fi tted LE profiles are shown in 

The bottom panels show the 
lit curves (left) and integrated 



the top panel of Figure 16 



corresponding effective Tig : 
spectrum (right). 

Note that while the LE profile fit is very good in all 
cases, the effective light curves are different, and thus 
the integrated spectra differ significantly in [O I] and 
[Ca II] line strengths. The reason for these differences 
is the degeneracy between the dust inclination and the 
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Fig. 14. — Effective light curves and integrated SN 1993J spec- 
tra varying the dust width, cr^ (top row), the dust inclination, a 
(middle row), and the seeing (bottom row) for LE2521. For a given 
row, only one parameter varies, with the other par ameters are set 
by the fit described above and shown in Figure |10| The black line 
indicates the original SN 1993J light curve and the corresponding 
integrated spectra. 
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dust width. For decreasing dust inclination (i.e., the 
dust inclination gets closer to the tangential of the LE 
paraboloid), the projected light curve is "stretched." 
This effect can be nullified by decreasing the dust width, 
resulting in a similar fit to the LE pr ofile with signifi- 
cantly different parameters (see Figure 15 1. 

For inclinations close to 90°, the LE profile is domi- 
nated by the dust profile and not by the projected light- 
curve profile. Such a situation manifests itself in a LE 
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Fig. 16. — Top: LE2521 profile and model fits for varying in- 
clination angles, a. Bottom: Effective light curves (left) and the 
corresponding integrated spectra (right). The black line indicates 
the original SN 1993J light curve and the corresponding integrated 
spectra. 

profile th at i s symmetrical (see green line in top panel 
of Figure 16 1 in contrast to the asymmetrical LE profile 



shape (see red line in top panel of Figure 16 1 that still 
reveals indications of the original light curve. Similarly, 
the fitted er^ is essentially constant for a ~ 90°. 

The tests in this section have shown that an incorrect 
dust inclination can systematically bias the interpreta- 
tion of SN properties from LE spectra. We discuss the 
sources of uncertainty in the determination of the incli- 
nation in Appendix IE] 

5. SN 1987A LIGHT ECHOES 

5.1. Observations 

Beginning in 2001, the SuperMACHO Project mi- 
crolensing survey employed the CTIO 4 m Blanco tele- 
scope with its 8K x 8K MOSAIC imager (plus its atmo- 
spheric dispersion corrector) to monitor the central por- 
tion of the LMC every other night for 5 seasons (Septem- 
ber through December) . The images were taken through 
our custom "VR" filter (A c = 625 nm, 5X = 220 nm; 
NOAO Code c6027). SN 1987A and its LEs are in field 
sm77, which was observed with exposure times between 
150 and 200 seconds. We have continued monitoring this 
field since the survey ended. The difference imaging re- 
duction process was identical to the KPNO Cas A imag- 
ing. 

5.2. Profile Modeling 

The ring-shaped LEs associated with SN 1987A have 
previously been studied in great deta il, mapping the 



dust st ructure surrounding the SN (e.g., Sugerman et al. 
|2005b[ ). Unlike historical SNe, a complete spectral and 
photometric history of the SN is available for SN 1987A. 
This historical database allows an observed LE spectrum 
to be modeled unambiguously against a spherically sym- 
metric model for the explosion. In addition, the high 
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Fig. 17. — Top: Apparent motion of a single position angle of the 
LE system of SN 1987A. The uniform motion of the LE on the sky 
indicates a unique inclination for the dust structure. The superlu- 
minal motion of the LE is 4.09 i 0.02 "/year. Bottom: The above 
(p, t) data is binned when observations are taken within 10 days 
of each other, and converted to the (z,p) plane using equation |Al| 
The resulting fit (solid line) gives an inclination of a = 30.4 ± 0.8"' 



signal-to-noise ratio of the SN 1987A LEs make them 
ideal for our tests. 

We consider an example LE from SN 1987A corre- 
sponding to a unique position angle on the sky. The spec- 
trum associated with this LE was take n using the Gemini 



Multi-Object Spectrograph (CMOS; |Hook et al.||2004 l 
on Gemini-South, using the R400 grating. The spectro- 
scopic slit was 1.0" wide and 2.16" long, with the slit 
tilted to be tangential to the observed LE ring. The up- 
per panel of Figure 17 shows the apparent motion of the 
LE monitored over more than a year. The motion on the 
sky is constant on such time scales (4.09 ± 0.02 "/year), 
indicating constant inclination of the dust structure. To 
determine the angle of inclination, a, we translate (p, t) 
into (z,p) as described in Appendix [E] resulting in a 
dust inclination of a = 30.4 ± 0.8° (see lower panel of 



Figure 17). 

A hign signal-to-noise ratio difference image is used 
to model the effect of the dust on the spectroscopic 
slit. Input parameters are the inclination (a = 30.4°), 
observed seeing for the image (1.22"), and delay time 
since explosio n (ct = 19.9 years). The photometric data 



of SN 1987A (|Hamuy et al. || 1988[ |Suntzeff et al.||1988 



Hamuy & Suntzctf 1990) are used to model the LL profile 



observed on the sky, with the best-fitting model shown in 
the left panel of Figure [18] The dust structure model is 
able to recreate the observed LE profile exceedingly well. 
The window function and resulting effective light curve 
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associated with the spectroscopic slit (convolved to the 
seeing at time of spectrosco py) are shown in the middle 
and right panels of Figure |18[ respectively. The aver- 
age dust width associated with the spectroscopic slit is 
a d = 4.26 ±0.06 ly (1.31 ±0.02 pc). We note that at the 
inclination and dust widths observed here, the model is 
less sensitive to uncertainties in a. Testing of the model 
shows that at such inclinations, uncertainties in a on the 
order of 5° have little impact on the windo w fun ction and 
resulting integrated spectrum (see Section 3.2). 

The impact of the scattering dust is clear when fitting 
the model to the observed spectrum. To integrate the 
spectral history of SN 1987A we take advantage of the 
extensive spectral database f rom SAAO and CTIO obser- 
vations after the explosion ( Mcnzics et al.||1987[ |Catch-| 
pole et al.||i~987l [19881 IWhitelock et al.||10l&l |Gatchpole| 
et al.|1989[|Whitelock et al.|1989| ?; 7)7 Figure [19] shows 
the observed LE spectrum (black), the model spectrum 
(red), and the full light-curve integrated spectrum (blue 
dashed). Note that the observed LE spectrum includes 
strong nebular emission lines that are not associated with 
the SN. Both the model spectrum and the full light-curve 
weighted integrated spectrum show a good fit for features 
blueward of 5500 A. However, the lower panel of Fig- 
shows that the spectrum redward of 5500 A, and 
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19 



in particular the Ha feature, are best fit by the spectrum 
corresponding to the best-fit LE profile model. From Fig- 
ure fl8 it is clear that the strength of the Ha feature in 



the integrated spectrum is highly dependent on the win- 
dow function. At late times, the spectra of SN 1987A are 
dominated by strong Ha emission. The excellent fit to 
the LE spectrum of SN 1987A validates our method of 
determining the effective light curve from the LE profile. 

5.3. Impact of Dust Parameters on Observed Spectrum 



Similar to the analysis performed in Section |4.3| we 
vary the three free parameters of our model to determine 
their impact on the observed LE spectrum for SN 1987A. 
Figure 20 shows the effective light curves (left column) 
and integrated SN 1987A spectra (right column). The 
top, middle, and bottom rows show the model varying 
the dust width (ad), dust inclination (a), and seeing, re- 
spectively, while holding the other parameters constant. 
To emphasize the differences in the integrated spectra, 
the figure focuses on the Ha feature. The advantage of 
the SN 1987A observations over those of Cas A is that 
there is a complete linking between the LE and the fi- 
nal integrated spectrum via observation, since we have 
a priori knowledge of the SN creating the observed LE 
profile. 

The dependence on the integrated spectra of SN 1987A 
are similar to that of SN 1993 J (see Figure 14 1, where 
increasing the dust width, inclination, or seeing all share 
the result of stretching the weight function closer to the 
wf = 1 case, resulting in an effective light curve ap- 
proaching that of the original SN 1987A light curve. Fig- 
ure [20 allows us to rank the impact of the dust parame- 



ters on the observed LE spectrum. Seeing has the least 
significance on the spectrum, although we note that poor 
image quality should still be taken into account. Dust 
inclination has a large effect, although the sensitivity is 
small if the dust is close to perpendicular to the observer. 
The dust width has the greatest impact on the spectrum, 



which is particularly noteworthy since this quantity has 
been ignored until now. For SN 1987A, the Ha peak 
height increases by about 20% if the dust thickness is 
10 ly instead of 5 ly, indicating a sensitivity that exceeds 
that of the other parameters. 

Figure [20} as well as the SN 1987A spectrum presented 
in Section~ |5.2| clearly show the danger in interpreting 
spectroscopic data from LEs without proper modeling. 
Any interpretation of LE spectra based on relative line 
strengths must be accompanied by a careful analysis of 
the dust properties associated with the LE. 

6. DISCUSSION 

In the previous sections, we have presented the LE 
model, the dependency of the model on various parame- 
ters, and have applied that model to multiple real-world 
cases. In this section, we discuss how we can use our new 
understanding of LEs and the scattering dust to perform 
additional observations that will provide new insight into 
historical SNe. 

6.1. Verification of Model Parameters And Systematic 

Biases 

The analysis of SN 1987A LEs has the added benefit 
that the validity of our dust models and their impact 
on the observed LE profiles and spectra can be tested, 
since the spectrophotometric evolution of the SN itself 
has been monitored e xtensively. W e show in Section [5] 
and in more detail in Sinnott et al. (2011) that the ob- 
served LE profile, apparent motion, and spectra can all 
be brought in excellent accordance with what is predicted 
using the SN 1987A spectrophotometric library and our 
dust model locally parameterized as a planar filament 
with a Gaussian density profile. This shows that at least 
in the case of SN 1987A, the Gaussian dust model is a 
good approximation of the true dust profile. 

In general, dust structures are filamentary in nature, 
which is easily visible in its effect on the shapes of LEs, 
which often show twists and furcations. However, on the 
few arcsecond scale, most of the shape and profiles of 
LEs change slowly over months and even years. For ex- 
ample, the dust width, o^, for LE3923 and LE2521 is 
nearly constant over 20" (see second panel from the bot- 
tom in Figure 12) and smoothly increasing, respectively. 



In order to account for small changes in the dust param- 
eters, we split the slit into 1" sections and fit the dust 
parameters separately for each of these sections. 

If there are significant deviations from the Gaussian 
dust density model (e.g., if there is substructure in the 
density on small scales) , the window function can poten- 
tially be biased, which in turn, can cause spectral differ- 
ences that would be interpreted as intrinsic to the SN. 
We can test and guard against these systematic biases: 
LEs often come in "groups"; i.e., there are often several 
distinct LEs separated by tens of arcseconds for LEs of 
Galactic SNe. The scattering dust filaments likely belong 
to the same dust structure, but are tens of light years 
apart. Therefore any substructure in these filaments is 
uncorrelated. For LEs of the same source event only tens 
or hundred of arcseconds apart, the physical properties 
of the source event are the same. An upper limit on 
the systematic bias can be set by the differences in the 
observed spectra which are not accounted for by the dif- 
ferences in fitted dust parameters. In all data obtained 
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Fig. 18. — Left: observed profile for the SN 1987A LE, with the best-fit model profile over-plotted. Middle: resulting window function 
associated with the 1.0" wide spectroscopic slit. Right: effective light curve from the model compared with the original light curve of 
SN 1987A. 
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Fig. 19.— Top: Observed LE spectrum of SN 1987A (black) with 
the best-fit model (using the window function from Figure [18| inte- 
grated spectrum (red) and the full light-curve weighted integrated 
spectrum (blue). Bottom: Detail of the spectra near the Ho fea- 
ture. Note that the narrow lines dominating the emission peak of 
the profile are strong nebular emission lines, and are not associ- 
ated with the SN. The best-fit model integrated spectrum fits the 
observed LE spectrum well in this region, while the full light-curve 
weighted integrated spectrum is a poor fit. 

by our team, we have not found any evidence that there 
are significant biases in our analysis due to substructure 
in the dust or deviation from the dust model. This is 
particularly important for the two novel methods that 
we describe in the next two sections. 

If LEs from two different SNe are reflected off the same 
dust, then the dust thickness measured from our method 
should be the same for both sets of LEs. Such a sce- 
nario is a direct test of the model fits. Fortunately, this 
situation occurs for LEs from Cas A and Tycho, and a 
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Fig. 20. — Effective light curves (left column) and integrated 
spectra (right column) for SN 1987A, varying the dust width, 
(top row), the dust inclination, a (middle row), and the seeing 
(bottom row). For a given row, only one parameter varies, with 
the other para meters are set by the fit described above and shown 
in Figure |18| 

comparison of the derived values will be presented in a 
future work. Since this situation was discovered in a rel- 
atively small and low filling factor survey, it is reasonable 
to expect other such regions of the sky exist where dust 
is at the intersection of two historical SN LE ellipsoids. 

6.2. Temporally Resolving Light Echo Spectra 

One of the most exciting opportunities LEs present is 
to use them to temporally resolve the SN spectrum: If 
a slit is aligned perpendicular to a LE profile instead of 
parallel, different parts of the slit probe different epochs 
of the SN light curve. 

As an example, we imagine a SN LE where the SN 
is a carbon copy of SN 1993J, the dust filament is at 
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z = 400 ly with an inclination of 0° and a width of 0.02 ly. 
With a Hubble Space Telescope-like (HST) PSF FWHM 
of 0.05", the projected light curve would have minimal 
distortion. If a spectrum is taken with the slit perpendic- 
ular to the LE arc (i.e., parallel to the line from the LE to 
the SNR) with 0.05" seeing, then the spatially resolved 
spectrum (i.e., moving along the axis perpendicular to 
the wavelength direction) corresponds to a t emp orally 
resolved spectrum. The left panel of Figure [2l] shows 
the effective light curves for this situation if spectra are 
extracted from 1" long parts of the slit. The offset from 
the peak of the LE are 1.0", 0.0", and -1.0" along the 
p-axis. The right panel shows the corresponding inte- 
grated spectra. The three spectra show clear differences, 
in particular the spectra with Ap ff sct = 1.0 is very blue 
since it is dominated by the first peak of the light curve. 

Aligning the slit perpendicular to the LE can signif- 
icantly reduce the amount of light within the slit, and 
therefore, this method is best applied to the brightest 
LEs. To maximize the temporal resolution of the LE 
spectra, the LE should have the following properties: 

1. The dust filament is very thin — on the order of 
0.01 ly. 

2. The inclination is close to the LE ellipsoid tan- 
gential. This "stretches" the projected light curve 
along the p-axis. 

3. The LE is bright. 

4. The PSF of the spectroscopic observation is small. 

Even though there are no known Galactic LEs that 
fulfill all of these criteria, we note that the existing sur- 
veys are far from complete. Additionally, even in non- 
ideal conditions, some temporal information could be ob- 
tained. The ability to temporally resolve the SN spec- 
tra will allow us to follow the evolution of the explo- 
sion, specifically measure velocity gradients of spectral 
features, and potentially probing the abundances of ele- 
ments in the ejecta. 

6.3. Constraining Supernova Light-Curve Shapes 
Under similar circumstances nee ded to temporally re- 



solve the LE spectra (see Section 6.2 1, we will also be 
able to spatially resolve the temporal characteristics of 
the source event light curve. The most important re- 
quirements are that the dust width has to be thin, the 
inclination close to the LE ellipsoid tangential, and the 
PSF needs to be significantly smaller than the width of 
the LE profile. Note that in contrast to temporally re- 
solving the LE spectra, constraining the LE shape can 
be done with LEs that have a significantly lower signal- 
to-noise ratio. 

Unfortunately, there is a degeneracy between the dust 
width and the light curve parameters (e.g., peak width 
and height; peak ratio if double-peaked) which can be 
broken if different areas of the LE have different dust 
widths. Since the light curve must be the same for all 
LEs in a given direction, any difference in the shape must 
be from the dust. If this degeneracy cannot be broken, 
LEs can still provide interesting constraints. First, the 
light curve can only be broadened by the dust, so any 
profile will provide some constraint on the light-curve 



shape. Additionally, by comparing the profile to known 
SN light curves, the properties of the SN can be further 
constrained. 

The light-curve shape of a SN can give important in- 
sights into the nature of the SN explosion. SNe la, for 
example, show a clear relationship between the width 
of their light curves and their intrinsic brightness (e.g., 
Phillips 1993 1. Another example is if the shock break-out 
of a core-collapse SN is resolved in the light curve, which 
would provide a measurement (or strict constraints) on 
the radius of the progenitor star. We will discuss this in 
more detail in the following paragraphs. 

After the gravitational collapse of the core of a mas- 
sive star, the core bounces, sending a shock wave out 
through the outer layers of the star, exploding the star 
( Woosley fc Weaver|1986 ). When the shock wave reaches 
the surface of the star, a flash (at optical wavelengths, 
it has a therma l spectrum) is generated that decays ove r 



hours to days (Colgate 19741 Klein fc Chevalier 



This has been seen fo r SNe UP (|Cezari et al. 
Schawinski et al 



2008), 



^2008 

but was best observed recently 



1978). 



with SIN lb 2008L) QSoderberg et al.||2008| |Modjaz et al. 
20091. As the stellar envelope cools, the optical emission 
lades away. Thi s subsequent fading was first detected in 
SN 1987A (e.g., |Woosley et al.|1987]|Hamuy et al.|1988[ ), 
which had a shock break-out cooling phase that was un- 
usually long. This cooling phase has been det ected for a 
handful of other SNe, notably SN 1993J (e.g 
et al.||1994h, a SN lib similar to Cas A. 



Richmond 



The duration, luminosity, and SED of the shock break- 
out phase depends on a handful of parameters s uch as the 
presence of a st ellar wind and the ejecta mass (Matzner 
fc McKee|1999|), but is most depen dent on the progenitor 
radius ( Calzavara & Matzner 2004 1 . The peak brightness 
of the subsequent fading is dependent on the ejecta mass, 
kinetic energy, and progenitor radius, while the timescale 
for the fading is proportional to both the radius of th e 



photosphere and its temperature (Waxman et al. 



2007) 



One caveat is that interaction of the shock with shells of 
material outside the star can produce prompt emission 
similar to a shock breakout. 

The disadvantage of using LEs to observe the shock 
break-out is that if the resolution of the image is poor or 
the dust thickness is too large, it is impossible to distin- 
guish between a shock break-out, subsequent fading, and 
normal radioactive-decay powered SN light curve, and a 
somewhat broader radioactive-decay powered SN light 
curve. For Cas A, where we might expect to see a shock 
break-out, we are unable to separate the light curve into 
multiple components in ground-based images with mod- 
erate seeing (< 1"). However, as shown in Figure [22| if 
the dust thickness is sufficiently small, we should detect 
a luminous shock break-out with the resolution of HST, 
if it exists. This will allow us to constrain the radius 
of the progenitor and its mass loss history, linking the 
physical properties of a star with its SNR. 

7. CONCLUSIONS 

In recent years, LEs of SNe have offered the rare oppor- 
tunity to spectroscopically classify the SN type hundreds 
of years afte r the light of the e xplosion first reached the 



Earth (e.g., Rest et al. 2008a). Recently, we have ob 



tained the first spectroscopic observations of a historical 
SN from different lines of sight using LEs (Rest et al. 
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Fig. 21. — Left: Effective fight curves for LE2116 with a 1" slit and offsets, Ap ff act , of —1.0", 0.0", and 1.0" shown in red, blue, and 
cyan, respectively. Right: Integrated spectra corresponding to the effective light curves in the left panel. 

the LE are is also a significant factor. The image PSF 
smears out the observed LE and thus also impacts the 
observed LE spectrum. Taking all of these factors into 
account, we found that the observed LE spectrum is not 
the light-curve weighted integration of the spectra at all 
epochs, but rather the integration of the spectra weighted 
by an effective light curve. This effective light curve is 
the original light curve convolved with a window func- 
tion, and it is different for every LE location. In general, 
the most significant difference between the effective light 
curves and the original light curve was that later epochs 
tend to be truncated in the effective light curves. 

We used fits to the LE profiles of Cas A and SN 1987A 
to test the consistency and veracity of our model. We fit 
the model to the observed LE profiles, where the three 
free parameters were the LE profile peak height, LE posi- 
tion, and - physically most important - the dust filament 
width. The model reproduced the observed LE profiles 
very well, and we found that it was possible to determine 
all factors influencing the observed LE spectrum a priori 
just from images, without using the observed spectrum. 

We constructed window functions and effective light 
curves for the Cas A LEs using SN 1993J as a template. 
We found a much better match of the SN 1993J spec- 
tra template to the observed Cas A LE spectrum if we 
used the effective light curve instead of the real light 
curve for weighting. The most significant difference was 
in the late-phase emission lines [O I] AA6300, 6363, [Ca II] 
AA7291, 7324, and the Ca II NIR triplet, consistent with 
our expectation fr om the model . Alth ough we do not 
have access to the Krause et al. (2008) images, we find 
that their Cas A spectrum is consistent with an inte- 




-1 

Ap in arcsec 

Fig. 22.— Simulated LE profiles of SN 1993J assuming differ- 
ent scattering dust widths. The scattering dust filament for the 
simulations is at z = 400 ly with an inclination of 0°, and the 
PSF FWHM is 0.05", similar to what can be expected from space- 
based imaging. The black dashed line is the projected SN 1993J 
light curve. 



20101. These observations provide the unique opportu- 



nity to examine a SN from different directions, with the 
light for each LE coming from slightly different hemi- 
spheres of the SN. In this paper, we have demonstrated 
that to properly separate differences in the observed LE 
spectra caused by intrinsic asymmetries in the SN from 
those caused by scattering dust, seeing, and spectro- 
scopic slit position, it is necessary to properly model the 
latter. 

Throughout this paper, we have shown how to both 
determine the properties of the scattering dust and their 
influence, as well as that of the PSF and spectroscopic 
slit position, on the observed LE spectra. We found that 
the dust width and inclination are the dominant factors, 
especially when the dust filament is thin. These two 
factors can be degenerate with respect to the LE pro- 
file shape; therefore, it is essential to determine the dust 
filament inclination independently using the LE proper 
motion from images taken at two (or more) epochs. The 
slit width and its misalignment with the major axis of 



grated spectrum weighted with a light curve truncated 
at 80 days. Our expectation is that the majority of dif- 
ferences between the Cas A LE spectrum and the full 
light-curve weighted SN 1993J spectrum are due to scat- 
tering dust properties. 

We also constructed a window function and effective 
light curve for a SN 1987A LE which constitutes a special 
test case since the light curve and spectral data of the 
SN are available. With this data, we constructed a LE 
spectrum from the original SN and compared it to the 
observed LE spectrum. We found that the constructed 
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LE spectrum correctly predicted the line strength of the 
Ha line, while a simple integration weighted by the light 
curve was a poor fit. 

We presented several additional investigations that can 
be attempted with LEs in the near future. Two of these 
studies, temporally resolved spectra and spatially re- 
solved light curves, require thin dust filaments. Until 
now, the assumption of thick scattering dust prevented 
such studies. 
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APPENDIX 
LIGHT ECHO PROFILES 



The LE equation ( |Couderc|p39| , 
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relates the depth coordinate, z, the LE-SN distance projected along the line-of-sight, to the LE distance p perpendicular 
to the line of sight, and the time t since the initial (undelayed) arrival of SN photons. Then the distance r from the 
scattering dust to the source event is r 2 — p 2 + z 2 , and p is 

p= (D-z)shi7, (A2) 

where D is the distance from the observer to the source event, and 7 is the angular separation between the source 
e vent and the sca ttering dust, which yields the three-dimensional position of the dust associated with the arclet. 



Tylenda (20041 elegantly derives the apparent expansion of a LE ring caused by a thin plane-parallel sheet of dust. 
In a very similar manner, we derive how the dust filament inclination influences the observed LE profile. Because of 
the symmetry around the z-axis, we define our coordinate system so that the dust sheet is normal in the (p, z) plane. 
This allows us to do all our calculations in the two-dimensional (p, z) space. We define the dust sheet as 

z = z d - ap, (A3) 

where a — tana. We define a as the angle of the dust sheet with respect to the p axis, where positive angles go 
from the positive p axis towards the negative z axis. The dust sheet intersects with the z axis at z = Zd- The two 
intersections of the dust sheet with the reflection paraboloid defined in Equation |A1| are then 



p±(t, Zd, a) = —act ± £, where (A4) 



£= y(l + a 2 ){cty+2z ct. (A5) 
Now consider a LE caused by a source event with a peak at to at an angular separation 7 from the source. Assu ming 



a di stan ce D to the source event, we can the n ca lculate the position of the scattering dust (po, Zq) using Equation Al 



and A2 The dust sheet defined in Equation A3 goes through (po,zo) if 

z d = z +ap . (A6) 

There are two special dust inclinations. We denote the inclination angle of the dust sheet tangential and orthogonal 
to the LE ellipse as ay and a±, respectively. Then the slopes of the dust sheets are an — tan ay, and a± = tanai, 
and we can calculate them as 

— v-(i+£)" 1/s -£ < A8 > 

The intersection of the dust sheet with the LE ellipsoid is then at 

fO for a = 90° 

Po = { P+ for a > ay . (A9) 
{p for a < ay 

Now consider a dust sheet with a finite width defined by a density D(Ar), and a light curve / SOU rcc(Ai), where Ar 
is the distance orthogonal to the dust sheet, and At = t — t (see Section IS] for a more detailed description of the dust 
model D(Ar)). The flux from the source event at time t is reflected by the part of the dust sheet at a distance Ar to 
the center of the dust sheet, and makes a relative contribution to the LE signal of 

/ LE = L»(Ar)/ Bource (Af) (A10) 

at the position 

_ (p+(h + At, z d + Ar/cosa,a) if p = p + (All) 
PLE \p-{to + At, z d + Ar/cosa, a) if p = p_ ' 

The reason for choosing + and — is based on the choice made for po- For all practical purposes, Ar and At are 
small, and therefore the sign should be the same. Exceptions are a « 0.0 and a as an. We define 

Ap = p LE - p . (A12) 

By examining the relative contribution for all values of Ar and At and using Equation |A10| - |"A12| we can calculate a 
cumulative grid of fluxes, F£ E (Ai, Ap). 
The model must now account for seeing. We describe the data PSF as a Gaussian with a FWHM of x arcsec and 



Gaussian width of ctpsf = x/(2y / 2 ln(2)) ps x/2.35. Then the cumulative grid of fluxes convolved with the seeing is 

E^E(At,A P ')^^exp(^ 



F LE (&t, Ap) = V Fl E (At, Ap') J^ exp ( {Ap n ^ P ' f ) . (A13) 
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Integrating over At returns the LE flux profile with respect to p: 

P LE (Ap) = J F LE (At,Ap)dt. (A14) 

DEPENDENCE OF THE LIGHT ECHO SPECTRUM ON LIGHT ECHO PROFILE 

We now examine the LE flux for a spectroscopic observation. Consider a slit positioned at p s with a effective 
width of er^cff- The reason that this is the effective slit width and its relation to the real, physical slit width a s is 
explained below in Appendix [C] We define the offset between the slit and the LE measured parallel to the p axis as 
Ap ff sct = p s — pq. The slit then subtends from 

Aftnin = Ap offsGt - <7 StC g/2 to (Bl) 

A/o max = Ap offset + a s , cS /2. (B2) 
We can then calculate the relative contribution, H(At), of the source event at time At to the spectra with 

H(At)= [ F LE (At,Ap)dAp. (B3) 

JA Pxair , 

H(At) is also the effective light curve. We define the window function w(At) as the fractional difference between the 
effective light curve and the real light curve, 

w(At) = f somcc (At)/H(At). (B4) 

The observed LE spectrum Sle(A) is then 

Sle(A) = J S(X,At)H(At)dAt, (B5) 

where S(X, At) is the spectrum of the source event at the given phase At. 

DEPENDENCE OF THE LIGHT ECHO SPECTRUM ON SLIT POSITION 

The position and inclination of the slit with respect to the LE has a major impact on the observed spectrum. We 
illustrate this in Figure [81 We define a coordinate system (p,u), where u is orthogonal to p through Ap = 0.0 in the 
plane of the sky. For a dust filament that is locally planar and perpendicular to the (p,z) plane (i.e., parallel to u), 
the LE is aligned with the u axis. However, if it is tilted with respect to u by an angle of Si, then the LE will also be 
tilted by the same angle (see gray shaded area in Figure [8]) . 

Now consider a slit (red box in Figure [8| that is misaligned with the LE by an angle 62, and also offset along the p 
axis by As. We split the slit into subslits, indicated with the dotted blue lines in Figure [8j For each of these subslits, 
we determine the LE profile in Figure [8j The offset Ap offsets of the subslit i centered at Ui to the LE is 

Ap s S et,i = ^ sin((5 2 ) + As, (CI) 

where 62 is the angle between the slit and the LE, and As is the offset of the center of the slit to the LE (see Figure [8]). 
The effective slit width cr S)e g depends on the physical slit width <j s and is 

c s , c ff = cr s /cos((5 2 ). (C2) 

The effective light curve Hi and window function uii for each subslit i can be calculated plugging Ap ff set) i into 



Equation Bl and B2 The total window function w B m for the full slit is then the weighted average of the subslit 
window functions, where the weight is the flux f w that falls into the subslit is 

/«,«= / " ma *P LE (Ap)dApwith (C3) 

JAp min 



-2 



w slit = — (C4) 

L^di J u 



> w . 1 



The integrated spectrum can then be calculated by using w s ut in Equations B4 and B5 

DUST MODEL 

For simplicity in our simulations, we choose a box-car dust density profile D(Ar). A more realistic density profile is 
a Gaussian, and therefore when fitting observed LE profiles we use a profile 
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where <Jd is the width of the Gaussian, and c is a normalization factor. Note that the physical dust density profile 
width cr^phys is different from the nominal fitted Ud if <5i 7^ 0. They are related through 

0"d, p h ys = a d cos(Si). (D2) 

DUST FILAMENT INCLINATION 

We refer the reader to a more thorough and detailed discussion about the apparent motion and its dependence on 
the dust filament inclination by Rest et al. (in prep.). However, we present an initial treatment of the process here. 
We define the dust filament inclination a as 

a = t & n-i(- Z 44^44) (El) 

with ti > ti, and with z calculated from Equation |A1| for a measured p and given t. Note that a dust filament in 
the plane of the sky has a — 0°, and dust filaments tilted away from the observer have positive a. If more than two 
epochs are available, fitting a straight line through (p, z) increases the accuracy of a. 
There are several sources that can introduce errors into the measurement of the dust inclination, specifically: 

• Poisson uncertainty in the apparent motion of the LE, especially for the fainter LEs 

• A small number of epochs 

• Intrinsic differences in the filament inclination with p 

The first two mainly increase the Poisson noise in the inclination. The last one, however, is the most significant one 
since it can introduce systematic uncertainties. If the dust filament has additional substructure, such as twisted dust 
on small scales, then the inclination of the scattering dust changes on small scales and a correct constraint on the dust 
inclination at the time of spectroscopy can only be done if there are several epochs close to the time of the spectroscopy. 



